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Abstract 

A pitchfork bifurcation of an (m — l)-dimensional invariant submani- 
fold of a dynamical system in R m is denned analogous to that in R. Suffi- 
cient conditions for such a bifurcation to occur are stated and existence of 
the bifurcated manifolds is proved under the stated hypotheses. For dis- 
crete dynamical systems, the existence of locally attracting manifolds M+ 
and M-, after the bifurcation has taken place is proved by constructing 
a diffeomorphism of the unstable manifold M. For continuous dynamical 
systems, the theorem is proved by transforming it to the discrete case. 
Techniques used for proving the theorem involve differential topology and 
analysis. The theorem is illustrated by means of a canonical example. 

1 Introduction 

Pitchfork bifurcations bear the name due to the fact that the bifurcation dia- 
gram for a one-parameter family in R looks like a pitchfork. Pitchfork bifurca- 
tion for a fixed point in R has been widely studied. In R, sufficient conditions 
for the occurrence of a pitchfork bifurcation of a non-hyperbolic fixed point are 
stated for instance in ^2 ^] . A generalization of the result in R is given by 
Sotomayor's theorem ^U] for a pitchfork bifurcation of a fixed point in R™. An- 
other generalization for pitchfork bifurcations is that for a periodic orbit |1(J|. 
Analytical discussions of pitchfork (or pitchfork type) bifurcations can be found 
for particular classes of dynamical systems, e.g.,|2], where a quasi-periodically 
forced map is studied. Interesting numerical analyses of pitchfork bifurcations 
can be found in Ej • 

An algorithm to compute invariant manifolds of equilibrium points and pe- 
riodic orbits is presented in [8]. It is important to study invariant manifolds 
in order to know the global dynamics of a system. The classical pitchfork bi- 
furcation concerns a fixed point (invariant codimension-1 submanifold) on the 
real line. From a mathematical viewpoint, it is therefore natural and impor- 
tant to investigate higher dimensional extensions of this theorem to invariant 
codimension-1 submanifolds of a Euclidean m-space. Accordingly we ask under 
what conditions an invariant manifold of a (discrete or continuous) dynamical 
system undergoes a pitchfork bifurcation? A fairly complete answer to this ques- 
tion is provided in this work. We give sufficient conditions for the occurrence 
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of a pitchfork bifurcation of a compact, boundaryless, codimension-1, invariant 
manifold in R m . We obtain readily verifiable criteria for identifying such bi- 
furcations, and illustrate the use of these criteria in an example. Techniques 
used for proving the theorem involve differential topology and analysis and are 
adapted from Hartman 4 , Hirsch et al. 5 and Shub |12| . 

2 Definitions 

Let M be a codimension-f , compact, connected, boundaryless submanifold of 
R m . By the Jordan-Brouwer separation theorem Pj, M divides R m \M into 
an outer unbounded region and an inner bounded region. We shall study C 1 
functions F : U x (—a, a) — > R m , where U is an open neighborhood of M in 
R m and (—a, a), a > 0, is an open symmetric interval of real numbers, ft shall 
be assumed in the sequel that each of the maps Fa : U — + R m , \fi\ < a, is a C 1 
diffeomorphism and that M is F M -invariant, i.e. F„(M) = M. 

Definition 2.1. With M and F^ as above, we say that F M is side-preserving if 
for every x in the inner bounded region, F^x) also lies in the inner region. 

Definition 2.2. With M and F^ as above, we say that F^ is side-reversing 
if for every x in the inner bounded region, F^x) lies in the outer unbounded 
region. 

Note that if F^ is a diffeomorphism of a neighborhood of M and leaves M 
invariant, then F^ is either side-preserving or side-reversing. Observe also that 
in the case that F M is side-reversing, it is not possible for U — R™. Analogous to 
the definition of a pitchfork bifurcation in K, we define a pitchfork bifurcation 
of invariant manifolds in R m as follows. 

Definition 2.3. Consider a discrete dynamical system in R m given by x n +i = 
F fl (x n ). Let M be an invariant manifold for all fi £ (—a, a). If < fio < a is such 
that M is locally attracting (repelling) for /i < 0, M is locally repelling (attract- 
ing) for \i > /io and in addition two locally attracting (repelling) -invariant 
diffeomorphic copies of M, viz., M_ and M + appear in a small neighborhood 
of M for /i > fiQ, then we say that M has undergone a pitchfork bifurcation at 
Mo- 

In the definition above, it does not matter what happens in the interval 
(0, /io). It is typically assumed that the interval (0, /io) is small. In R, coincides 
with /io, since the manifold under consideration is just a single point. But 
for higher dimensions, not all points on the invariant manifold may undergo a 
change in stability at the same value of the parameter \x. When /i > /io, all 
the points have changed stability and two new invariant, diffeomorphic copies 
of the original manifold (of opposite stability) appear. 
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3 Pitchfork bifurcation theorem for discrete dy- 
namical system 

We consider (one-parameter families of) discrete dynamical systems given by 

x n+1 = F(x n ,fx) (1) 

where, x n £ R m for every n £ N and [i £ (—a, a) C R. Additional properties 
of F(-,[i), also denoted as -F M (-) are described later. Let M be a compact, 
connected, boundaryless, codimension-1, C 1 submanifold of R™, which is F^- 
invariant V/i G (—a, a). Denote a tubular neighborhood of M as N(a) — {x £ 
R m : <i(x, M) < a, a > 0}, where d is the standard Euclidean distance function. 
Assume that a is sufficiently small so that the e-neighborhood theorem [3] can 
be applied and N(a) C U. This means that every element x £ A^(a) can be 
uniquely represented as a; = (r, y) where j/ = tt(x) £ Af is the point on M closest 
to x and r £ [— a, a] is the signed distance in the outward normal direction 
between x and M. We also assume that F )Ji (N(a)) C N(a). This enables 
us to write F^ in component form as F^ = (/^,<7p) where / M : N(a) — > R, 
,9;t = 7T -F/i : AT(a) — > M, and /^(ic) is the signed distance from ^(x) to 
F[t(x). Observe that F~ 1 (i f1 p (iV(a))) = -/V(a), so that F^ 1 can be written in 

component form as F^ 1 = (Z^,^). 

We shall use standard notation for derivatives and partial derivatives of 
functions. For example, the derivative of F^ : N(a) — > R m will be denoted by 
DF^, and represented as the usual m x m Jacobian matrix 



DF ll {r,y) = 



D r U(r,y) D y f^(r,y) 
D r g^(r,y) D y g^(r,y) 



where the entries are submatrices representing the partial derivatives such as 



df 

D rU(r,y) = -g^{r,y) and Dyg^r^y) 



dyj 



(m— 1) X (m— 1) 



We use | • | for the Euclidean norm of an element of a Euclidean space or 
the associated norm of a linear mapping (matrix) between Euclidean spaces. 
The symbol || • || denotes the supremum norm of a function taking values in 
a Euclidean space or in a space of linear transformations of Euclidean spaces 
taken over an appropriate set, which is sometimes indicated as a subscript of 
the norm. 

If the rate of change (with respect to r) in the normal component in 
the radial direction r is strictly less than 1 in absolute value, the manifold M 
will be locally attracting. This is stated mathematically in statement (ii) of 
Theorem 13. II which follows. Similarly to have M locally repelling, we require 
that \Drfnl be greater than one as in statement {Hi). Statements (iv) and (v) 
describe locally attracting properties in a neighborhood away from M, which 
is where our new bifurcated manifolds M_ and M + will reside. Properties 
(vi) and (vii) are obtained analytically and are needed in order to establish 
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the existence of manifolds M_ and M + as graphs of a fixed point (Lipschitz 
function) in a Banach space. The last hypothesis (statement (viii)) provides 
boundedness and equicontinuity properties that enable us to bootstrap Lipschitz 
homeomorphisms of M with M + and M_ up to C 1 diffeomorphisms. These ideas 
shall become clear as the proof unfolds and after the remarks following the proof. 

Theorem 3.1. With and M as above, suppose that the following statements 
hold. 

1. is side-preserving for every fi 6 (—a, a). 

2. sup \D r fp(r,y)\ = \\D r f p \\ N ( a ) <l for every iie{-a,Q). 

(r,y)eN(a) 

3. 3 < /i* < a such that inf \D r fu(0, y)\ > 1 V/i G (//*,a). 

4-3 < ot\ < a such that ||D r / /J (r, v)\\a < 1 V/i G [0, a), where A = 
{x G W l : a\ < d(x, M) < a}. 

5. 3 x '■ [0j a ) ~ > ^ continuous with < x(/i) < ai and K(fi) := {a; G 
K m : X(M) < d(x,M) = d(x,y) < a} such that F^(K(fj,)) C JsT(/i) V/i G 
(fJ,*,a). Furthermore c(/i) := \\D r f^\\ K ^ < 1 V/i G {^*,a). 

6. c*(/i) := + ||-Dr<yk(p)) + H-Dk/^IU^) < 1 for each fi G 
(/i*,a), where \\.\\k(ij.) is defined to be the sup norm over K([i). Here 
(//ijff/j) denotes the inverse map F~ . 

7. (H-Dr/^lk^) + ll^/^lk( A1 ))(ll- D r5plk(p) + ll-D^lkcrf) ^ 1 /° r eac/l M G 

(/^.AO- 
'S- O-(M) : = ll^r/ M |k( M )(2||I)r3Ml|K(M) + ll I:) !/3Mlk( M )) + ||-D y /Ml|x( M )||A-ff M ||K( M ) < 1 

/or a/Z /i G (/i*, a). 

TTien for each \x G (/i*, a), 3 codimension-1 submanifolds M + (p) and M-(p) 
in K(p) such that both M + (p) and M_(/i) are F ^-invariant, locally attracting 
and C 1 diffeomorphic to M. M is locally repelling, and r > for all x = (r, y) G 
M+, and r < for all x = (r, y) G M_. 

Proof. Recall that AT(a) = {x = (r,y) G JR m : |r| = d(x,M) < a , y = 
7r(a;)}. Now F^(r,y) = (f^(r,y),g^(r,y)) in component form where, : 
N(a) — > R is the signed distance between i^(r, y) and g^(r, y) and gy, : N(a) — > 
M is the projection 7r o F fJi (y) of F^{r,y) on M. We shall break the proof up 
into a number of steps (claims). 

Claim 1. M is locally attracting for [i G (— a,0). 

Proof of CZaimQJ- Consider a point (ro,?/o) S iV(a). Let (r n ,y n ) be the point 
obtained by applying the n-fold composition of F^ with itself to (ro,yo). Then 

(r n ,y n ) = F M (r n -i,y n -i) = (.fii( r n-l>yn-l),9iJi(rn-l,Vn-l)) 
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implies that 



d((r n ,y n ),M) = d((r n ,y n ),Tr(r n ,y n )) = / M (r n _ X) y n -i). 
As M is F^-invariant, it follows that /^(O, y n -i) = for all n G N. So, 

\r n \ = \f^(r n ~i,yn-i)\ = \fn(r n -x, y n -i) - U(0,y n -i)\ = \ d f^ r ' Vn ~^ ||r n _i| 



by the mean value theorem. Thus |r„| < c|r„_i| < c"|ro|, where 



c = sup I d ^'( r ' v '> I < 1 by property (m). Therefore, r„ — > as n 



(r,y)&N{a) 

Consequently d((r n ,y n ), M) — ► 0. That is, for any initial point (ro,J/o) hi the 
neighborhood AT(a) of M, F"(ro,yo) converges to M. It follows that M is 
locally attracting for all /i G (—a, 0). 

Claim 2. M is locally repelling for fj, G (fj,+ , a). 

Proof of Claim 0' Following the same steps as above, we find that |r n | > 
c l r n-i| > l r n-i| whenever |r„| is sufficiently small owing to statement (Hi). 
Accordingly the iterates {x n } must eventually leave any sufficiently thin tubular 
neighborhood of M for /i G (/!*, a), which means that M is locally repelling. 

We now fix a fi G (/i* , a) and suppress /x in the notation for simplicity. To 
begin with, we shall prove the existence of M+ as an i^-invariant manifold 
homeomorphic to M. It suffices to prove the existence of M+, as the existence 
of M- can be established in the same way. Observe that M+ is invariant iff 
F(M + ) — M + . We shall seek M + in the form of the graph of a continuous 
function over M defined as 

M + =T il = {(iP(y),y):yeM}, 

where M+ C K = K(p), i\) : M -> M and ?/%) > for all y G M. Then for all 
y £ M, we have that (i>{y),y) € M + iff F(ip(y),y) = (ip(z),z) G M +) which is 
equivalent to 

/(^^jW.^WW.^^ (2) 

f is a diffeomorphism, hence F (ip(z), z) = {ip(y),y) which implies that 

(f(1>(z),z),g(<<i>(z),z)) = (i>(y),y). (3) 

where F^ 1 = (/,<?). Combining equations @ and ©, we find that M + is 
invariant iff if; satisfies the functional equation 

^(z) = fU(g('iP(z),z)),g(ij(z),z)\ (4) 
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Let Lip(A, B) denote the set of all Lipschitz functions from A to B. Let £(ip) 
denote the Lipschitz constant of a Lipschitz function ip, and L^ = {(ip(y),y) ■ 
y G M} denote the graph of ip. Now define the set 

X :={ip£ Lip(M,R+ U {0}) : C(tp) < 1,1^ C K}. 

Claim 3. {X, \\ ■ \\k} is a Banach space. 

Proof of Claim\$ Let {ip n } be a Cauchy sequence in X. Then for all n G N 

we have that %P n : M M+ U {0}, £(>„) < 1 and G Here L^ 

denotes the graph of ip n over M. Since the sequence {ip n } is Lipschitz, every 
ip n is continuous. Now M is compact and M+ U {0} is closed, so the set of 
all continuous functions from M to K + U {0} with sup norm forms a Banach 
space. Moreover, if ip n — > ip as n — > oo, it is clear that ?/> is also Lipschitz, with 
Lipschitz constant not greater than one. Since K is closed, K contains all its 
limit points. Therefore G K for all n implies that 

Urn T^ n = Urn (ip n (y), y) = (ip(y),y) = I> G K, 

n — >oo n — >oo 

hence X is a Banach space. In view of we define an operator T on X as 
follows. 

H^){y) :=fU(g^(v),y)),g^(y),y)\ (5) 

Claim 4. C X. 

Proof of Claim^ Let z = g(ip(y),y). Then f(ip(z),z) — signed distance be- 
tween F(ip(z),z) and n(F(ip(z),z)). If ^(z) > 0, then f(ip(z),z) > since F 
is side-preserving. So .^(VO is indeed a function from M to ]R + U {0}. ^(ip) 
is continuous since it is a composition of continuous functions. Now it follows 
from the mean value theorem and the definition of X that 

m){Vi)-H^){yi)\ = \fmzi),z 1 )-f(i;(z 2 ),z 2 )\ 

< W^WkMzi) - ip(z 2 )\ + \\D y f\\ K \z x - z 2 \ 

< (\\%\\K + \\D y f\\K^ \Z!-Z2\ 

where \\.\\k — sup{\.\ : (r, y) G K}. The above inequality follows because ip is a 
Lipschitz function with Lipschitz constant < 1. Also 

\zi-z 2 \ = \g{ip(yi),yi) -g(tp(y 2 ), y 2 )\ 

< HA-slk Vl>{vi) - ^(2/2)1 + ll-Dj/fflk \yi-V2\ 

< [ \\D r g\\K + \\D v g\\ K ) \y1-y2l 
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The two inequalities obtained above, together with property (vii) imply that 



\H4>){yi) - H4>){y2)\ < 



< 



dr 

x IS/i — 2/1 1 

\yi — 2/2 1 - 



K + \\D v f\\ K )[\\D r g\\ K + \\D y g\\ K 



Therefore, F(%j)) E Lip(M, K+ U {0}) and C(T(ip)) < 1. We will now prove that 
T is a contraction mapping. Using statement (vi) and the mean value theorem, 
we compute that 

\H^){y) - Hi>2){y)\ 

= \f(M(Mv),v)),9(Mv),v)) - f(*h(g{Mv),v)),9{Mv),v))\ 

df 



< 



< 



1^1(5(^1(2/), y)) - i>2(g(i>2(y),y))\ 



K 



dr 

+\\Dyf\\ K \9(My),y)-9(My),y)\ 

df 



dr 



\i>i{g{i>i{y),y)) - i>i(g(i>2(y),y))\ 



K 



df 



dr 



\i>\ {g{ip2 (y),y)) - ^2(5(^2(2/), y))\ 



K 



< 



< 



■\\Dyf\\ K \g{Mv),v) - 9iMv),v)\ 

df 



(g(i>i(y),y) -g{^2{y),y) + \\i>i - V>2| 



K 



dr 

+\\Dy.f\\K\g(My),y) - g(My),y)\ 

df 



dr 
df 



(WDrgUU, - V2II + IIV-i + H^/lkll^lkllV'i - M 



K 



dr 



K 



(l + \\D r g\\ K ) + \\D y f\\ K }Ui-M\- 



This is true for all x G K . Hence it is true for the supremum with x taken over 
K; therefore, we obtain the relation 



ll-FfrM-JWII < 

11^0-^2)11 < 



cJVi-^II 



\D r g\\ K ) + \\D y f\\ K 



IIV'l II 



where c* = || ||k(1 + ||A-s||.k-) + IIA//Ik is such that < c* < 1 by hypothesis 
(vi). This implies that Tty) e Lip(M,R+ U {0}) and that C(T(%1>)) < 1. We 
note here that the invariance of M implies that f(0,y) = f(0,y) = for all 
(0,2/) £ M. Accordingly D y f(0,y) = D y f(0,y) = whenever a; = (0,2/) £ M, 
which means that both H-D^/Hic and can be made as small as we 

like by choosing a sufficiently thin tubular neighborhood of M. Note that if 
-P(^(y)>2/) = (ip( z ), z )> tnen (^(j/))2/) — I p ~ 1 (' t l J ( z )j z )i an d it follows that 



^(y) = f(i>(z),z) 



and 



2/ = ffW'W, z). 
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Now consider ip € X. By definition, we have 



I> W = {(f(i>)(z),z):z£M} 

= \ \f('l>(§Mz),z)),g(il>(z),z)),z\ :zeM 



This implies that 

= {(fMy),y),z) :zeM} 

= {ifMv),v),g(iKy),v)) : gWy),y) e M} 

= {(/(M»).ffWl/).V)):»eM} 

= {F(ip(y),y) : y £ M}. 

We know that {^{y),y) <E if for all y e M and F(iiC) C if. This implies that 
r^(^) £ if, thereby proving the claim that ^(X) C X. Hence J 7 : X — > X is a 
contraction mapping with respect to the sup norm on X . 

Since T is a contraction on a complete metric space X, it has a unique fixed 
point in X owing to Banach's fixed point theorem. Let <fi be the fixed point 
of T. Then </> G Lip(M,R+ U {0}) with Lipschitz constant £(0) < f, and 
satisfies the functional equation (@J. Therefore, 

ct>{z)=fL{g{<t>{z),z)),g{<t>{z),z)\ (6) 

Claim 5. M+ exists and is locally attracting. 

Proof of Claim We now define M + as the graph of as follows 

M + =r* = {(0(y),y);yeM}, 

where is as above. This proves the existence of M + . That M + is locally 
attracting follows directly from its definition as the graph of a fixed point (func- 
tion) of a contraction mapping. 

Claim 6. M+ is homeomorphic to M . 

Proof of Claim: Let H : M -> M+ be defined as H(y) := (0(y),y). Then 
ii is injective, surjective and continuous. ii _1 exists and is also injective and 
surjective (bijective). Since M + is compact, H~ x is also continuous. Hence the 
manifold M + is homeomorphic to M. 

Claim 7. The function is a class C 1 map. 

Proof of Claim: We know that is the solution to the functional equation 

®, hence 0(z) = G Lip(M,R+ U {0}), and 

£(0) < 1. We will inductively construct a sequence of C 1 functions tp n which 
converges to 0. Then using the Arzela-Ascoli theorem, we will prove that is 
C 1 . The details are as follows. 
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Choose ipi to be a positive constant such that C K. By construction, 
ipi is C 1 and C{\j)\) = 0. Now suppose ip n is defined and that ip n is C 1 with 
£{i>n) < 1- We define ip n +i inductively as, 

i> n+1 (z) = H*Pn)(z) (7) 
= / (^Pn(g(lpn(z), Z)), g(lp n {z), z)\ 

Let h n : M — > ]R™ 1 denote the function (ip n ,Id), where Id denotes the identity 
map on the second coordinate. That is, h n (z) = (ip n (z),z). Then h n is C 1 by 
the induction hypothesis and the fact that it is the composition of C 1 maps, and 

ipn+i(z) = f ° h n o g o h n (z). 

Here we have used the fact that both / and g are C 1 since F and F -1 are C 1 
diffeomorphisms . 

The sequence of functions {tp n {z)} converges uniformly to <j)(z) since <f> sat- 
isfies the contractive functional equation Q- The Jacobian of ip n +i evaluated 
at z is the 1 x (m — 1) matrix or the gradient vector of ip n +i given as 

Dip n+1 {z) - I3/^ n (5(/i n (^))^£»/i„(g(/i„(z)))£>g(/i n (z))-D/i n (^), (8) 

owing to the chain rule. Moreover, 

ljj n +l(z) = T(4) n )(z) 

by construction. Hence, C(ip n +i) < 1. Since, V'n+i is differentiable, this implies 
that ||D-0 n +i(z)|| < 1 for all z G M. By induction, {Dijj n (z)} is a sequence of 
continuous functions, uniformly bounded by 1. 

We will now prove the equicontinuity of {Dip n (z)}. The techniques used 
below are actually global versions of the methods employed by Hartman 0] 
for local invariant manifolds, and the role of the Lipschitz property follows an 
approach used by Hirsch et al. |Hj and Shub to study hyperbolic invariant 
manifolds. 

For any function /?, we define A/3(z) := (3(z + Az) — (3(z). When Az is such 
that Az < minjj, ^ D gjr^p; w^ }i we wm show that || ADip n (z)\\ < t(8) for 
all n, where r depends only on 5 and is such that r(S) — > as 5 — > 0. The 
desired result will be proved by induction as follows. 

For any i5 > 0, we define quantities rj(S) and t(S) as 

V (S) = sup | ||D r /(r + Ar, y + Ay) - D r f(r, y)\\, 

\\D y f(r + Ar lV + Ay) - D y f(r, y)\\, \\D r f(r + Ar,y + Ay) - D r f(r,y)\\, 
\\D y f(r + Ar,y + Ay) - D y f(r, y)\\, \\D r g(r + Ar,y + Ay) - D r g(r, y)\\, 
\\D y g(r + Ar,y + Ay) - D y g(r, y)\\, \\D r g(r + Ar,y + Ay) - D r g(r,y)\\, 

\\D y g(r + Ar,y + Ay) - D y g(r,y)\\ : (r, y) G N(a), \ Ar\, \Ay\ < s\ 
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and 



t{6) = 2(\\D r f\\ K + \\D y f\\ K + \\D r g\\ K + \\D y g\\ K ) 

1 — a 

where a < 1 is as defined in property (viii). It is observed that r)(6) converges 
to as S approaches 0. 

Recalling that tpi is defined to be a constant, we have ADip\{z) = 0, and this 
implies that ||ADr/>i(z)|| < t{5) for all Az. Suppose that \\ADip n (z)\\ < t(S) is 
satisfied whenever Az < min j<5, nr , -n — fn n -n }. Now, 

— L ' \\ D rg\\K + \\D y g\\ K 1 ' 

Dip n+1 (z) = DxCxBxA 

where 



D = [ D r f(ip n (g(ip n (z),z)),g(ip n (z),z)) D v f(ip n (g(tp n (z) 7 z)) 7 g(tp n (z),z)) 
is a 1 x m matrix, 



C 



D^n{9{lpn{z),z)) 
Im-1 



mx (m— 1) 



and 



B= [ Dg{^ n {z),z) 



DM*) 

Im-l 



I (m-l)xm ' 



mx(m— 1) 

Multiplying the four matrices above and taking into account the block matrix 
notation, it follows that Dip n+ i(z) can be expressed in the following simpler 
form. 



Dip n +i(z) 

= D r f(ip n (g(ip n (z), z)),g(ip n (z), z))Dtp n (g(ip n (z), z))D r g(ip n (z), z)Dip n (z) 
+D y f{ip n (g(ipn(z),z)),g(il) n (z), z))D r g(ip n (z), z)Dip n (z) 
+D r f(i>n(g(ipn(z), z)), g(ipn(z), z))Dtp n (g(ip n (z) , z))D y g{ip n (z) , z) 
+D y f(tp n (g(tp n (z),z)),g(ip n (z),z))D y g(tp n (z),z). 

Each of the four terms added above is a 1 x (m — 1) vector. We will now 
estimate the quantity j ADip n+ \(z)\\. Using the definitions, we find after a 
straightforward calculation that AD\p n+ i(z) = Dtp n+ i(z + Az) — Dtp n+ i(z) 
can be written in the form 
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AD^ n+1 (z) 

= ^D r f(^ n (g(^ n (z + Az),z + Azj), g($ n (z + Az), z + Az)) 

xDxP n (g{iP n (z + Az), z + Az))D r g(^ n (z + Az), z + Az) 
xDip n {z + Az) 
-D r f(ip n (g(ip n (z), zj), g(ip n (z), z)) 

xDip n (g(ip n (z), z))D r g(tp n (z), z)Dip n (z)^ 

+ ^D v f(i/) n (g(i/) n (z + Az), z + Azj), g(^ n (z + Az),z + Az)) 

xD r g(iP n (z + Az), z + Az)DtP n (z + Az) 
-Dyf(tp n (g(t(; n (z),z)),g(t(; n (z),z))D r g(tf; n (z),z)Dtp n (z)^ 

+ ^D r f(^ n (g(^ n (z + Az), z + Az)),g(ip n (z + Az), z + Az)) 

xDip n (g(ip n (z + Az), z + Az))D y g(i) n {z + Az), z + Az) 
-D r f(ip n (g(tp n (z), z)), g(ip n (z), z))Dip n (g(ip n (z), z))D y g(tp n (z), z) j 

+ ^D y f(4> n (g(ip n {z + Az), z + Az)),g(ip n (z + Az), z + Az)) 

xD y g(tp n (z + Az), z + Az) 
-D y f(ip n (g(ip n (z),z)),g(ij n (z),z))D y g(il) n (z),z)y 

We denote the four bracketed terms above as T\, T 2 , T 3 and T4, respectively. 
For instance, 

T 2 = D y f{^ n (g{ib n (z + Az),z + Az)),g{ib n {z + Az),z + Az)) 
xD r g(ip n (z + Az), z + Az)Dtp n (z + Az) 
-Dyf{ip n (g{ip n {z), zj), g(ip n {z), z))D r g(ip n (z), z)Dip n (z). 
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Adding and subtracting appropriate terms yields, 
T 2 

= I.D y f(iJ n (g(ip n {z + Az), z + Az)), g{^ n {z + Az), z + Az)) 

xD r g(ij) n (z + Az), z + Az)Dtp n (z + Az) 
-D y f(ip n (g(ip n (z), z)),g(ip n (z), z))D r g{ip n {z + Az), z + Az) 

xDtp n (z + Az) 



+ i Dyf(ip n (g(ip n (z), z)), g(ip n {z), z))D r g(ip n (z + Az), z + Az) 
xDij) n (z + Az) 

-D y f(ip n (g(-il) n (z), z)), g(ip n (z), z))D r g(t/j n (z), z)Dip n {z + Az] 



+ 1 D yf{i } n{9{4>n{z),z)),g{^ n {z),z))D r g{'ip n {z),z)Dij n {z + Az) 
-D y f(tp n (g(ip n (z), z)), g(tp n (z), z))D r g(ip n (z), z)Di/j n (z) 



Using the triangle inequality and the above definition of the quantity ?y(<5), we 
obtain 

||T 2 || < V (S)\\D r g\\ K \\Dij n \\ K + \\D y f\\ KV (6)\\Di; n \\ K 
+ \\D y f\\ K \\D r g\\ K \\ADi; n \\, 

where the first term is valid only when Az < „ n ~„ — f M n - M — and Az < 5. 

J — \\D r g\\ K + \\ D v9\\K — 

(The estimate on Az is obtained by applying the chain rule and the mean value 
theorem to the first term in T 2 .) Since £(ip n ) < 1, it follows that ||Z?i/>„|| < 1- 
By the induction hypothesis, ||AZ?-0„|| < t(5). This implies that 

||T 3 || < v(6) (\\D r g\\ K + ||U w /||jc) + \\D y f\\ K \\D r g\\ K r(S). (9) 

Using similar analyses, we obtain 

\\Ti\\ < r)(6)(\\D r g\\ K + \\D r f\\ K ) + 2\\D r f\\ K \\D r g\\ K r(6), (10) 

||T 8 || < vm\\D r f\\ K + \\D y g\\ K ) + \\D r f\\ K \\D y g\\ K r(5) (11) 

and 

\\T4<vm\D y g\\K + \\D y f\\ K ). (12) 
Combining equations © - (|12fl gives, 

||A^n+i|| < 2{\\D r g\\ K + \\D y g\\ K + \\$£\\ K + \\D y f\\ K } V (6) 

+ {\\§\\K{2\\D r g\\ K + \\D y g\\ K ) + \\D v f\\ K \\D r g\\ K r{5)) . 
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Substituting the definition of a in the above inequality, we find that 
IIAU^ill <t(S)(1-o) + ot(S), 

which proves that || ADip n+ i \\ < t(S) whenever Az < min{<5, ^ D g \\ K + \\ D ^ K }■ 
Thus, we have proved by induction that ||AD^„(z)|| < r(S) (whenever Az is 
sufficiently small) for all n. The quantity r(<5) is such that, t(8) — > uniformly 
as 5 approaches 0. Hence, the sequence of functions {D?p n (z)} is equicontinu- 
ous. Since the sequence {Dtp n (z)} is a uniformly bounded and equicontinuous 
sequence of functions on a compact set M, it follows from the Arzela-Ascoli the- 
orem that there exists a subsequence D?p nk (z) which is uniformly convergent on 
M. Let p(z) be the uniform limit of D?p nk (z) as k — > oo. Since we know that 
ip n converges to cf>, this implies that p — D(f>. That is, <p is differentiable. Also, 
since Dip n (z) is continuous for every n and the convergence is uniform, we find 
that p is also continuous. That is, D<p(z) is continuous. This implies that <fi is 
class C 1 . 

Hence, the map H : M —* M + defined earlier as H(y) :— (4>{y),y) is a C 1 
diffeomorphism. Thus we have proved that the manifold M + is diffcomorphic 
to M. 

Analogously, one can prove that Af_ is diffeomorphic to M . This proves that 
M has undergone a pitchfork bifurcation at into a pair of locally attracting 
invariant manifolds M+ and M_, each diffeomorphic to M , for each p € a). 
Thus the proof is complete. □ 

There is also a side-reversing version of Theorem 13.11 that can be proved in 
a completely analogous manner; namely 

Theorem 3.2. Let and M satisfy all the hypotheses of Theorem Vd '. 1\ except 
with being side-reversing. Then for each p £ (p+,a), there exist manifolds 
M-ip) and M+(p), both C 1 diffeomorphic to M, such that F^{M + ) = Af_ 7 
F^(Af_) = M+ and M_(/i) U M + (p) is F ^-invariant and locally attracting . 

In certain cases, the estimates in properties (vi)-(viii) of Theorems 13 . II and 
13. 21 can be combined into a single statement, as in the following result. 

Corollary 3.3. Let the hypotheses of Theorems \3.1\ and VS.IA be as above, except 
that properties (vi)-(viii) are replaced by the single estimate 

(ix) Pr/ M k||^k + (||^/p||A- + P a / M k)(^l + || J Dr^||A^ < 1 for each 
p G (/i* , a) . Then the conclusions of the theorem still follow. 

Proof. One need only observe that (vi)-(viii) follows directly from (ix). □ 

Remark 3.1. If the function F^ and invariant manifold M are of class C 2 , prop- 
erty (viii) of the above theorem is not essential, for in this case the equiconti- 
nuity of the sequence {Dip n } in the above proof follows from the mean value 
theorem. Of course, if one wishes to prove the existence of bifurcated C 2 diffeo- 
morphs of M, an analog of property (viii) involving second derivatives would be 
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necessary. Such an estimate, although rather complicated, can be obtained in a 
straightforward manner, and we leave this to the reader. If both the map and 
invariant submanifold are C k , with k > 2, it is not difficult to obtain a fcth order 
derivative analog of (viii) that would guarantee the existence of bifurcated C k 
diffeomorphs of M. 

Corollary 3.4. Let the hypotheses be the same as in Theorems Vd '. 1\ a,nd \3.tA with 

the following additional modifications: Property (v) is replaced by 

(if') 3x '■ [0, a) such that < xit 1 ) < a i> an d ^(Kdi)) C K{jjL), 

where K(fi) is as in (v) for every (i* < \i < a, and the following assumption is 

added. 

(x) For every fi G (/i*,a), f„(r,y) > r(< r) for (r,y) G (0,x(m)] x M and 
Fp(r,y) < r(> r) for (r,y) G [— %(/i),0) x M in the side-preserving (side- 
reversing) case. 

Then, in addition to the conclusions of Theorem \3.1\ and \'3. A we have the fol- 
lowing dynamical properties: The submanifold M+(/i) attracts all points x = 
(r,y) G (0, a] x M, and M_(/x) attracts all points x = (r,y) G [a,0) x M in 
the side-preserving case; and in the side-reversing case, N(a)\M is contained 
in the basin of attraction of M + (/i) U M_ (/x) . 

Proof. We shall verify only the additional result for M + (p) in the side-reversing 
case, since the proofs of all of the other cases are similar and require only 
obvious modifications. For the case at hand, it obviously suffices to show that 
the iterates of a point (ro, yo) with < ro < xO-O eventually wind up in K(p). 
Setting (r n ,y n ) — F™(ro,yo), it follows from (x) that {r n } is an increasing 
sequence of real numbers, which must exceed x(p) f° r n sufficiently large. Thus 
the proof is complete. □ 

Remark 3.2. It is natural to ask about the bifurcation phenomena that may 
occur when < fi < fi* and M has regions where \D r f^\ > 1 and regions 
where \D r f^\ < 1. If F M leaves all points of M fixed, one can readily prove the 
existence of "blistered" diffeomorphs of M using one-dimensional theory. The 
"blister" regions, where \D r fn\ > 1, have a pair of locally attracting copies of 
M manifested as inner or outer blisters on M, while the portion of M inside 
the blister is locally repelling. However, when merely leaves M invariant 
without fixing all the points, the situation is much more complicated and needs 
further investigation. 

Remark 3.3. A particularly useful feature of our main results, Theorem 13.11 
Theorem 13.21 ( and Theorem 15 . 21 appearing in the sequel), is that they are con- 
structive. The desired bifurcated manifolds can be determined to any desired 
accuracy by successive approximation. For example, to approximate M + (fi) in 
the side-preserving case, one simply starts with ipx as a positive constant so 
that its graph is in K(fi), and then computes successive approximations using 
the functional equation Q. The iterate ip n for n sufficiently large yields an 
approximation M n that can be chosen to be arbitrarily C 1 close to M + (/n), and 
the error can be estimated from the definition of the iterates. 
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4 Illustration of the (discrete) pitchfork bifur- 
cation theorem 



In this section, we illustrate Theorem 13.11 proved in Section [3] with a canonical 
example. Let A £ SO„(R), the special orthogonal group of real n x n matrices, 
comprised of orthogonal matrices with determinant 1. Define a linear map 
L A : K" -> K" as 

La(x) = ^4x. 

The map is an analytic (linear) diffeomorphism. Every (n — l)-sphere S a 
of radius a > is Lyi-invariant. That is, La{So) = S a where S a = {x 6 K™ : 
|x| = a} and a > 0. Note that Si denotes a sphere of radius 1, in the space on 
which La acts. If for instance, A <E SC>2(R) then La ■ R 2 — > R 2 and Si is the 
same as S 1 C K 2 . If ^4 e S0 3 (R) then La : R 3 -> K 3 and Si is the same as 
S 2 cl 3 . The subscript denotes the radius of the sphere and the dimension of 
the sphere is one less than the ambient space. 

Now define : [0, oo) — > [0, oo) to be a C°° function such that tr M satisfies 
the following properties. 

1. <7^ = in a small neighborhood of 0. 

2. ct m (s) > 1 for < s < f. 

3. o- M (s) = 1 - (s - l) 3 + m(* - 1) for I < s < |. 

4. cr M (s) < 1 for I < s. 

5. (so>(s))' = scr^(s) + cr^(a) > for fj, £ ±}. 

We fix a matrix A in S0„(R) and define F p : R" -> R™ as follows: 

F M (x) = a- M (|a;|)iA(a;) = o>(|x|)Ax. 

It is easy to see that is a diffeomorphism, and that it leaves Si invariant. 
That is, F /1 (Si) = S\. The discrete dynamical system governed by F M is 

x n+ i = Fn(x n ). (13) 

In the notation of Section |3 M = Si. Due to the symmetry of the sphere Si, 
every point in R™\{0} can be uniquely described as being a radial projection on 
Si, so the neighborhood N(a) is not restricted by the e-neighborhood theorem. 
However, due to the nature of tr M , we let a — -g and consider the neighborhood 
iV(|) = {x e M™ : x e [|, |]}. We now check that all the hypotheses stated in 
Theorem 13. II are satisfied. 

1. Observe that F^ is side-preserving for fi 6 ^g] since A preserves ori- 
entation and cr M is positive-valued. 
For this example, 

x 

r = x — 1 and y = - — -. 
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This implies that after a change of variables, F^{x) = ( | x \)Ax becomes 

F»(r,y) =a li (r + l)\x\A£ [ , 
F»(r t y) =a lt (r + l)(r + l)Ay. 

The property that A preserves length is used in obtaining the above ex- 
pression for F^, and again in finding / M and below: 

U(r, y) = |(T M (r + l)(r + l)Ay\ - 1 = (r + l)a„(r + 1), 

9ii{r,y) = Ay. 

This implies that 

df 

^- = (r + iy fl (r+l) + a IM (r + l), 
D y f(r, y) = 0, D r g(r, y) = and D tf5 (r, y) = A. 

2. sup |^| < 1 for all /i. € [-jg-,0) since the maximum 1 is attained at 
[X = as shown in Figure 2] 




Figure 1 : rvs^ for the canonical ex- 
ample for r G [—0.2, 0.2] as \x increases 
from ^ through 0. 

25 ° 



3. For this example, /i* = and inf | g° | > 1 for all \i S (0, ^]. The 
infimum is attained at /i = as illustrated in Figure 

4. For this case, o,\ can be chosen to be 0.15. As illustrated in Figure |3 
sup 1 1£| < 1, where A = {(r,y) : 0.15 < r < 0.2}. 

A 

5. K{n) can be chosen to be A and property (v) follows from property (iv). 

6. Properties (vi), (vii) and (viii) also follow from statement (iv) since 
D rgii (r,y) = 0, D„/(r,y) = and ||Z> w ff(r, y)|| = 1. 
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Figure 2: Plot of jj, vs ^£ for the 
canonical example for r — and in 
the interval [0, 




Figure 3: Plot of r vs ^ for the canon- 
ical example for r G [—0.2,0.2]. 



Theorem 13. II implies that S\ undergoes a pitchfork bifurcation at /z* = 0. This 
is indeed the case and for fi 6 (0, 1/25]: F u has three invariant spheres Sj_ 
Si and S'l+yj: where Si is locally repelling, and S±^rp and Si+^i are locally 
attracting. This is illustrated in Figures 0] El and 

Remark 4.1. The above example can be easily modified to illustrate Theorem 
EHl Define a map = R o F^, where R : R m \{0} -» R m \{0} is a smooth 
map such that R(x) — ^jjjp on the neighborhood N(^) of Si. Then G M is 
side- reversing, with all other properties the same as those of Fa- In this case, 
S± undergoes a pitchfork bifurcation at /i* = and for /i £ (0, jg] the invariant 
manifolds are Si and Si_/^ U Si +v /^-. Note that G M (Si_y^) = Si+^tz and 
G At (S'i + yjr) = S'l-yjr. 
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Figure 4: Any trajectory outside Si Figure 5: Any trajectory inside Si 
converges to Si. converges to Si. 

For /i < that is, before the bifurcation Si is locally attracting. Figures 01 and 
[5] are generated for fi — ^ in the canonical example. 




Figure 6: Any trajectory outside con- Figure 7: Any trajectory inside con- 
verges to Si+^p. verges to Sx-^a- 
For [i > 0, Si is locally repelling. Figures El and are generated for [i = 
Several points are iterated many times to generate these pictures. 



5 Pitchfork bifurcation theorem for continuous 
dynamical system 

In this section, we state and prove a pitchfork bifurcation theorem for continuous 
dynamical systems that is analogous to the result for the discrete case given by 
Theorem l3.ll The idea of the proof is to use the flow generated by a continuous 
system to reduce the problem to the discrete system covered by Theorem 13.11 
Consider a continuous dynamical system given by 

x = X(x,(j,) (14) 
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where x £ K m and [i E (—a, a) C R. The (1-parameter) vector field X(x,fi), 
also denoted as X fl (x), is assumed to be of class C 1 in its domain. 

Let <j>{t,x,n), which we also denote by </^(x), be the unique solution (flow) 
starting at x when t — 0, and M be a compact, connected, boundaryless, 
codimension-1, ^-invariant manifold of R m , which means 0*(M) = M for all 

< a. As in Sectional we define 

iV(a) = {x E M m : d(x, M) < a} 

as a tubular neighborhood around M, where the e-neighborhood theorem [2] is 
applicable. Any point x in the region N(a) can be written as 

x = (r, y) 

where r is the signed distance between x and the manifold M and y is the 
unique point of M closest to x. Recall that r is positive if x lies in the outer 
unbounded region of M. m \M and r is negative if x lies in the inner bounded 
region of W n \M. Again as in Section|21 the notion of inner and outer regions is 
obtained as an application of the Jordan-Brouwer separation theorem 3 . Note 
that r — when x lies on M. 

We shall assume that the vector field points into N(a) on dN(a) for every 
ii in the set (—a, a), which means that positive semi-orbits of (jl4|l that begin 
in N(a) can never exit this tubular neighborhood. Note that X^ = (R^^Y^) 
in r,y-component form. Analogous to Section |3J it follows that <^ maps N(a) 
into itself for all (t, fi) E [0, oc) x (—a, a). Now we can write the flow in terms 
of r and y components as 

<K*> 0, J/), M) = G°(*> ( r > ?/))/*)> ^(*) ( r > y)> 

where </>(t,(r,y), is the signed distance between <f>(t,(r,y), fx) and M, and 
^(t, (r, y), /i) is the normal projection of <f)(t, (r, y),fJ,) onto M. 

In order to obtain the estimates necessary to reduce the continuous case 
to the discrete case, we shall need to consider the derivative of the flow with 
respect to the initial condition x = (r, y). By D(j)(t, x, //), we mean the Jacobian 
matrix of <f> with respect to x defined as 



D<f>(t,x,fi) 



D r p{t,x,ii) D y p(t,x,Li) 
D r ip(t,x,/ji) Dyif>(t,x,fi) 



It is well known that the matrix Dcj>(t, x, li) is the unique solution of the initial 
value problem (see e.g. in Hartman 0]) 



<S> = D x X^(t,x,ii))§ = 



D r R^ DyR^ 



$ (15) 



$(0) 



where I m is the m x m identity matrix. 

We shall make use of the following version of Gronwall's inequality. 
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Lemma 5.1. Consider the linear matrix initial value problem 



$ = r(i)$ 
*(0) - I m , 



(16) 



where $ = {4>ij), r = (7^) and T is a continuous matrix function on the real 
line M. Let cj>(t), <5>i{t), $//(*), $>ui(t), j(t), r/(t), T H (t) and T ln (t) be the 
submatrices defined as follows: 

0(i) = 11 (i), 7 (i) = 7 1 i(f), 

$/(*) = [012W,--- ,^lm(*)],rj(t) - [7l2(t),- •■ , 7 lm(4)], 



= [M*).-- ,0mi(*)f .r//(*) = [72lW 

022(4) ■•• 02m (t) 



*///(*) 



,T ln (t) = 



( Pm2 (4) ■■■ 0mm (4) 

where the superscript T denotes transpose, so that 



$j 
$11 $111 



7(4) r 7 (i) 
r w (t) r //7 (t) 



722(4) 

7m2(4) 



$J 



72m (4) 
7mm (4) 



0(0) = 1, $/(0) = 0, $//(0) = 0, */(0) - I m - 
numbers satisfying 

a, v, a 2 ,v 2 < s/4, 
and suppose that for some positive t*, 



Let a, v and s be positive 
(17) 



7(4) < -2s, |r z (t)|, \T n (t)\ <a, \T m (t)\<u, (18) 
whenever \t\ < t*. Then for all \t\ < t* the solution of \lb}) satisfies the estimates 



X-t 



(j(A + + 2s) 1 



x+t 



X-t 



10(4)1 <E {t) :=Kie 
\$ll(t)\ <E 2 (t) := a{\ - v)- 1 Kl e 

|*/(t)| < E x {t) := K ie A -* - a(X+2s)K 2 e 
\*in(t)\ <E 3 {t) := aiX-v)- 1 ^-* - K 2 e X+t , 



n 2 e x +\ 

-,_ X + t 



(19) 



where 



A+ := - 



(2s - v) 



1± 4 1 



4cr 2 



(2s-i/) 2 



(20) 



/si := 1 - <t 2 {(A+ + 2s)[(A+ + 2s) + a 2 (A_ - v)]}' 1 , 
k 2 := -a{(A_^)(A+ + 2s)[(A+ + 2s) + a 2 (A_ - v)]}' 1 , 
«i := a(A_^)[(A+ + 2s)(A_ - 1/) + cr 2 ]" 1 , 
75 2 :={l + ( T 2 [(A + +2s)(A_- i ,)]- 1 }- 1 . 
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Proof. It follows from equations H16fl - i|18|) m the hypotheses that 

\<P(t)\ < u(t), \$i(t)\ < v(t), \$n(t)\ < w(t), \$iu(t)\ < z(t) 

for all \t\ < t- kl where u, v, w and z are the entries of the 2x2 matrix initial 
value problem 



it 


V 




' -2s 


a 




u 


V 


w 


z 




a 


V 




w 


z 



(21) 



" "(0) v(0) 1 _ 
w(0) z(0) \ ~ 2 

The eigenvalues - one negative and denoted by A_ , and one positive and denoted 
by A + - of the constant matrix in (|21|l are easily computed and found to be given 
by (|20|l . Now l|21|l can be solved by elementary means to yield 

u(t) = E (t), v(t) = E^t), w(t) = E 2 (t), z{t) = E 3 (i), 

where Eo, E%, E 2 and E$ are as defined in l|19|) . Accordingly, we have verified 
the desired estimates, thereby completing the proof. □ 

Theorem 5.2. Let the vector field X : M. m x (—a, a) be C , and let M be a 

compact, connected, codimensions-1 invariant manifold for \1J$ for every /i € 
( — (j, (7 ) . Suppose that the following properties hold: 

1. X^ points into N(a) for all (x,/!) £ dN(a) x (—a, a). 

2. D r R(x,(i) < for all x — (r,y) in the neighborhood N(a) for all fi S 
(-a,0). 

3. There exists < /i* < a such that D r R(x, /i) > for all (x, /x) € M X 
(ft*, a). 

4- For each n £ (^*,a) there exists < ai(/i) < a and an s > such that 
X^ points into 

A(/i) = {x e E m : aiO«) < d(x,M) < a} 
on its boundary and D r R((r,y), [i) < —2s for (r,y) G A((i). 
5. Let a and v be positive constants such that a , v, a 2 , v 2 < s/4 
\\D y Rj A(M) , \\DrYj < a, \\D y Yj A(M) < v 
for all n S (/!*, a), and a, v are sufficiently small with respect to s so that 
E (t)(l + E 2 (-<))+ E^t) < 1, (22) 

(E (t) + £?i(t))(^(-t) + E 3 (-t)) < 1, (23) 
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E (t)(2E 2 (-t) + E 3 (-t)) + E 1 {t)E 2 (-t)) < 1, (24) 

where E , E\, E 2 and E 3 are as in Lemma \5.1\ for all fi G (/i*,a) and 
each 1 < t < 2. Then the invariant submanifold M is locally attracting for 
fi G (— a,0), and locally repelling for \i G (p+,a). Furthermore, for each 
fi G (/i*, a) £/iere exist a pair o/ C 1 diffeomorphs M + (p) and M_(/i) o/ 
Af in smc/i £/ia£ feot/i and M„(/i) are invariant for and 

locally attracting. 

Proof. Using the relation (r, y) = (R((r,y), /j,),Y((r,y), //)), it follows from the 
mean value theorem that 

r = i?((r, y). M ) = R((r, y).fj,) - R((0, y),fj,) 
= D r R((r„,y),^)r, 

where relies between and r, and we have used the property i?((0,y), [i) = : 
which follows from the invariance of M . Consequently, property (ii) implies that 
r < when r > and r > when r < 0, which means that trajectories tend 
toward M as t increases. Hence, M is locally attracting for each —a < fi < 0. 
Similarly, one can also use the mean value theorem to show that it follows from 
property (Hi) that M is locally repelling for fi G (/i*, a). 

From here on, we fix // G (/i* , a) and suppress it in order to simplify the 
notation. We shall first show that for each t G [1,2], the map T* defined as 

T*(a;) := 00, x), 

where ^ is the flow generated by the differential equation (|14|l . satisfies the 
hypotheses of Theorem lH.il As $ = D(f> satisfies the initial value problem (|15|l . 
Lemma 15. II is an ideal instrument for proving the desired result. 

Observe that from the form of the estimates i|19|) of Lemma 15.11 that for 
a given s > it is indeed possible to select positive numbers cr, v sufficiently 
small for estimates (I22H - I|24|) of property (iv) to obtain for all 1 < t < 2. This is 
with the understanding that we may assume without loss of generality that we 
are in T 2 (A), so that we can still take advantage of the initial norm estimates 
for the terms with arguments —t in l|22|) - (|24|l . which correspond to the inverse 
of T* owing to the group property of the flow <f>. Accordingly it follows from 
Lemma |5.1l and Theorem 13. II that for each /i G (/i*,a) and t G [1,2], the map 
T* has a unique pair of contractive invariant manifolds M± in A, which are 
C 1 -diffeomorphic with M . 

We shall now show that the manifolds M|j_ and M*_ are, in fact, the same 
for all t G K, and they are invariant for the entire flow <f>. It is enough to verify 
this for Ml, as the proof for M l _ is identical. Consider any rational number of 
the form q = 1 + — lying between 1 and 2. Then, by definition 

T q {M q + )=M q + . 

Applying the map T q , n times to this equation yields 

[T q ] n {M q + ) = Ml, 
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which by the additivity property of the flow becomes 

r (n+m)( M 9) =M\. 

But the unique contractive manifold for T 1 is Mi, and r(™+ m ) is an (n+m)-fold 
composite of T 1 with itself. Hence, 

T (n+m)( M l) = 

so it follows from uniqueness that 

Ml = Ml, (25) 

and this must hold for all rational numbers 1 < q < 2. 

It now follows from l|25|) . the completeness of real numbers, and the conti- 
nuity of the flow that Mj = Mj. for all t € [1,2], and 

T*(M[) = M[ (26) 

for all 1 < t < 2. Any t > 2 can be written as t = m + r, where m is a positive 
integer and r € [1,2]. Consequently, 

T f (Ml ) = T (m+T ) {Ml) = T m o T T (Ml ) 
= T m (Ml) = Ml, 

owing to H26|) and the definition of M^. Whence l|26() holds for all t > 1. In 
fact, it holds for all \t\ > 1 since 

T~ t (T t (Ml)) = Ml = T- l {Ml) 

whenever t > 1. 

Finally, for any e > 0, 

r £ (Afi) = rf^or 1 ^) = T< 1+e) (M+) = Mj. 

Thus T*(M{) = M| for all t > 0, and it therefore follows as above that the 
same is true for all i < as well. We denote the unique invariant manifold 
by M+. This yields the desired result that T\M + ) = M + for all tel, which 
completes the proof. □ 

Remark 5.1. For continuous dynamical systems, only the side-preserving case 
can occur. 

Just as in the case of a discrete dynamical system, we can obtain a more 
complete description of the dynamical systems in N(a) by making a minor 
additional assumption. 

Corollary 5.3. In addition to the hypotheses of Theorem 15. ill suppose that 
R({r,y), n) is positive (negative) for < r < a\(fi) (-ai(n) < r < 0) whenever 
(i 6 (/it*, a). Then M + (/i) attracts all points of N(a) with r > and Af_(/i) 
attracts all points of N(a) with r < 0. 

Proof. The additional property guarantees that the positive semi-orbits in N(a) 
not lying in M, eventually enter A(n) and are then attracted to M + (fi) or 
M-(fi). This completes the proof. □ 



23 



6 Conclusions 



We have proved that codimension-1, compact invariant manifolds in discrete dy- 
namical systems, undergo pitchfork bifurcations when the system satisfies suit- 
able conditions. The hypotheses of the theorem are easily verifiable estimates 
on the norms of partial derivatives of the function determining the discrete dy- 
namical system, which makes this result well suited to a variety of applications. 
When the bifurcation parameter /i is between and /x*, some portions of M 
may be locally repelling and some locally attracting (in the normal direction), 
so the proof of our theorem would need to be modified to handle this case, which 
is an interesting subject for future investigation. 

The case when the whole manifold M bifurcates into M_ and M + as /x 
increases through zero, corresponds to //* = 0. The fact that can be greater 
than allows for M to eventually bifurcate and does not impose the restriction 
that M bifurcate all at once. The theorem is slightly weaker than the theorem in 
one-dimension since the theorem does not completely determine the dynamics 
of the system in the region between a neighborhood of M and the neighborhood 
A of M- and M+. 

The pitchfork bifurcation in K is assumed to be one stable fixed point bifur- 
cating into two stable fixed points separated by an unstable fixed point. We have 
generalized this result to a compact, connected, boundaryless, codimension-1, 
locally attracting invariant submanifold of E m that becomes locally repelling 
and bifurcates into two locally attracting diffeomorphic copies of itself separated 
by the locally repelling manifold. The techniques we have used here should en- 
able us to obtain new results on higher dimensional versions of other types of 
bifurcations such as Hopf and saddle- node Hopf bifurcations (see e.g. 0). We 
plan to investigate these and related generalizations in the future. 
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